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Microscopic membrane elasticity and interactions among 
membrane inclusions: interplay between the shape, dilation, tilt 
and tilt-difference modes 
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Abstract. A phenomenological Landau elasticity for the shape, dilation, and lipid-tilt of bilayer membranes 
is developed. The shape mode couples with the sum of the monolayers' tilt, while the dilation mode couples 
with the difference of the monolayers' tilts. Interactions among membrane inclusions within regular arrays 
are discussed. Inclusions modifying the membrane thickness and/or inducing a tilt-difference due to their 
convex or concave shape yield a dilation-induced attraction and a tilt-difference-induced repulsion. The 
resulting interaction can stabilize 2D crystal phases, with the possible coexistence of different lattice 
spacings when the dilation-tilt-difference coupling is large. Inclusions favoring crystals are those with 
either a long-convex or a short-concave hydrophobic core. Inclusions inducing a local membrane curvature 
due to their conical shape repel one another. At short inclusions separations, a tilt comparable with the 
inclusion's cone angle develops: it relaxes the membrane curvature and reduces the repulsion. At large 
separations the tilt vanishes, whatever the value of the shape-tilt coupling. 

PACS. 34.20-b - 82.65Dp - 87.22Bt 



1 Introduction 

Bilayer membranes are formed by the self-assembly of am- 
phiphilic molecules in water of brine The aliphatic 
chains of the constituent molecules condense into an oily 
sheet that is shielded from contact with water by the po- 
lar heads of the molecules. Membranes can form various 
phases, ej* lamellar (L a ), vesicular (£4), or sponge (L3) 
phases 0|| . Being self-assembled systems with a con- 
served area, it is essentially the competition between their 
curvature energy and their entropy that determines their 
large scale behavior. In the standard macroscopic theory, 
membranes are modeled as structureless surfaces with a 
curvature elasticity [§,1). This description, which has the 
advantage to involve only two material constants, accounts 
for a large number of universal properties and behaviors 
of amphiphilic membranes |^,^,^]. 

On the other hand, a number of attempts have been 
made toward a more microscopic description of membranes 
|@,|[||,|l^,|ll|,^2|. The goal being to take into account vari- 
ous structural parameters of the bilayer, such as its thick- 
ness, the ordering of the chain segments, etc. As far as 
large scale properties are concerned, these extra degrees 
of freedom are irrelevant since they relax over microscopic 
lengths. Nevertheless, they can dictate important physical 
properties, such as adhesion behaviors, the short-range in- 
teraction between membrane inclusions, their aggregation 
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properties, phase behaviors, etc. The aim of this work is to 
construct an elastic model of membranes that connects mi- 
croscopic and macroscopic descriptions. Beside the stan- 
dard shape and dilation variables we shall consider as elas- 
tic variables the tilts of the lipids in both monolayers [|l3| 
|l4| . The model will be used to investigate the role of the 
monolayers tilts in the interaction between membrane in- 
clusions. 



2 Elastic model 

To construct an elastic model, one considers a distorsion 
free energy depending on a particular set of structural pa- 
rameters. Implicitly, this free energy results from the in- 
tegration over all the microscopic states compatible with 
these (fixed) parameters of the Boltzmann weight asso- 
ciated with the microscopic Hamiltonian of the system. 
In practice, based on the symmetries of the system, one 
writes an expansion in powers of the structural parameters 
and their gradients. The choice of the relevant parameters 
depends on which deformations can be imposed externally 
on the system. We shall consider four structural parame- 
ters: (1) the membrane thickness, which can be modified, 
e.g., by the presence of an integral protein, (2) the mem- 
brane average shape, in order to connect with the large- 
scale theory and because it can be excited by a conically 
shaped inclusion, and (3) the tilts of the molecules within 
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Fig. 1. a) Average membrane shape, h, and dilation, u, variables, b) Construction of the membrane average tilt m = |(p 
p (2) ) + Vh and tilt-difference fn = |(p (1) + p (2) ) variables. 



each monolayer, that can be independently excited by an 
inclusion with, e.g., a diamond-like shape. 

To simplify, we assume that the membrane undergoes 
only small deviations with respect to its flat ground state. 
We denote by h^\x,y) and h( 2 >(x,y) the vertical dis- 
placements (along z) of the chain-water interfaces of the 
upper and lower monolayer, respectively, with respect their 
positions in the flat unperturbed state (Fig. For fur- 
ther use, let us define the average shape h(x,y) and the 
membrane dilation by 



(i) 

(2) 



To construct the tilt variables, we introduce the vectors 
p^ 1 ) (x, y) and (x, y) defined in both monolayers as the 
projections onto the (x, y) plane of the unit vectors par- 
allel to the molecular direction and oriented from chain 
to polar head (Fig. |l|b). The tilts relative to the mem- 
brane normal are measured by m' 1 ' = p' 1 ' + V/i and 
V/i. Let us define the average tilt m(x,y) 



m 



(2) = n (2) 



and the tilt- difference m.(x, y) by 

m (l) _ m (2) 



m (l) + m (2) 



(3) 
(4) 



2.1 Shape and dilation distortion energy 



We start by constructing the most general quadratic free 
energy expansion in powers of hP~> and ft/ 2 ) and its first 
and second gradient, 



» 



(5) 



Here, a — 1,2 is the monolayer label, and the comma 
denotes partial derivation with respect to the coordinates 
x and y. We write the free energy as F = F^'+F^'+F^- 12 ' 
with all the interaction terms coupling h^ 1 ' and in 
_F(i2)_ The symmetry of the bilayer imposes invariance 
with respect to the transformation: 

h^^-h^ 2 \ /iP) _» . ( 6 ) 
Therefore, the most general quadratic form for is 

+athWh$ +a n hf h { f+a 6 h^ 2 

(7) 



-a 7 n tij ti tij 



summation over repeated indices being understood. The 
interaction energy, containing all the bilinear scalars, has 
the form 



F 



12) 



-b 3 hy^ +bi h$m+hh$h%. ( 8) 



Expressing now F in terms of h and u, we obtain the 
decoupled form F = Fh + F u , with 

F h = d\ h 2 + d 2 h t i h t i + d 3 h h,u + d± h 2 u 

+d 5 h ti j h^j (9) 

F u — e\ u + e 2 u 2 + e 3 u,j u ti + e 4 u t u + e 5 u u y u 

+e e u 2 u + e 7 Uij u y ij . (10) 

The new coefficients are related to the former by an in- 
vertible linear transformation. In this expression, several 
terms obviously vanish and other can be discarded: e\ = 0, 
since the minimum energy corresponds to u = by con- 
struction; g?i = d% = since F must be invariant under a 
translation. We shall set do = 0, since the tension of mem- 
branes usually vanishes |fl5[ [| . There is no reason however 
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to discard the term e3(Vw) 2 , which represents the energy 
density associated with a gradient of the membrane thick- 
ness. The latter term involves not only the extra coast 
of lengthening the chain-water interfaces but also that of 
modulating the stretching of the molecules chains. Con- 
siders a planar membrane with a thickness modulation 
at some wavevector q: its elastic energy may be well de- 
scribed by the term ex u 2 as long as qa <C 1 (a monolayer 
thickness), however the term cx (Vu) 2 should not be ne- 
glected when qa w 1. From this point of view the present 
model differs from those of Refs. J^,[n^[ll],[l2| that neglect 
the coefficient 63 ]lq| . 

We can now rewrite F in a more traditional way. Rela- 
beling the nonzero coefficients, and making use of hjj hjj = 
(V 2 /i) 2 — 2Det(/i,y), we arrive at 



Fu 



F„ 



k(V 2 /j) 2 + &Det(h tij ) , 

B u 2 + -A (Vw) 2 + a V 2 u + a' u W 2 u 



+ -k' (V 2 m) 2 + k' Det(u (j j) . 



(11) 



(12) 



Fh is simply the Hclfrich energy in which V 2 h is 

twice the mean curvature of the average membrane shape, 
and Det(/i i i J ) is its Gaussian curvature. The thickness 
variations, which are completely decoupled from the aver- 
age membrane shape, are described by an energy F u simi- 
lar to that of Refs. pQ , p| ,p^, however with two important 
differences: (1) there is a non vanishing term cx (Vm) 2 at 
lowest-order, (2) the bending constants k' and k' are dif- 
ferent from the Helfrich constants appearing in F^. 

To further simplify our model, we shall discard the 
terms proportional to a and a', since they can be trans- 
formed to boundary terms by integration by parts; we 
shall also discard the terms proportional to k' and k' , in 
order to keep only the leading-order saturation terms. We 
are left (at the moment) with 



and the exchange of and m/ 2 ). Writing all the linear 
and quadratic scalars yields 



G (q) = Aim 

G (12) 

r<(a) 



) I a (a) (a) . , (a)2 

+ A 2 m\ m\ '+A 3 mi^ 

a (a) (a) a (a) (a) 

+A4 m\ J m] J + A 5 m\ J m) A ' 
B\ m) 



B3 m 



C\ m 



1) (2) 



T3 (1) (2) 

B 2 rn ■ fn ■ 

z m 3,3 



,(!)_(2) 



h3 
a) 



1,3 



D (1) (2) 

B 4 m\j m) > 

(«), 



u.i + C2 m) i' u + C3 m 



(15) 
(16) 
(17) 



(-l) a Di m[ a) h ti + {-l) a D 2 m\f h 
+{-i) a D 3 m^ h dj + (-lr^mg fc,«(18) 



As previously, several terms can be discarded: A\ — C 2 — 
D2 = 0, since we assume no spontaneous splay of the 
tilt; D\ = since the minimum energy is still achieved 
with zero tilts when the membrane is rotated. We can 
also discard the terms with coefficients A5 and B4: inte- 
grating them by parts merely yields boundary terms (and 
a renormalization of A3 and B 2 ). 

In terms of the variables m and m, we can write the 
total tilt energy as G = G m + G m h + G — h G^ , with 



l h3 



1 2 
G m = -t mi rrn + fci m i t + k 2 m it j 

G m h — <^3 m i,i h,jj + <^4 m i,j n Aj 

G m = -t' fhi rhi + k'x m 2 ^ + k' 2 rh iti m it j 
G~ u = cfhi u,i + ci m iti u tjj + c 2 m hj u t ij 



(19) 
(20) 
(21) 
(22) 



We can finally discard the terms with coefficients c 2 and 
c?4 by integrating by parts, and neglect the term with co- 
efficient ci as a higher-order coupling term. 



F= {V 2 K) 2 + RT>et(h,ij) + ^B u 2 + ^A (Vu) 2 (13) 



2.2 Tilt distortion energy and coupling terms 



2.3 Total distortion energy 

In vectorial notations and after some simple manipula- 
tions, the total distortion energy, i.e., F + G, can be writ- 
ten as Hhm + H — , with 

1 li iv urn ' 



We expand the distortion energy associated with the tilts 
of the molecular orientation in powers of 



m 



The tilt gradient Tn^j is a non-symmetric second rank and 

tensor. We write the tilt free energy as G = G^ + G^ + 
G^ 12 ) + Gint where all the terms coupling and m^ 2 ) 
are in G^ 12 - 1 and all the interaction terms coupling the 
tilts and the membrane shape or dilation are in Gmt- 
The interaction can be divided into four contributions: 



m 



(«). 

i,3 ' 



(14) 



int — 



(1) 



"11 



G 



(1) 



( 2 \ 

G h , each term containing 



all the contributions bilinear in either rrS^ or to' 2 ' and 
u or h. Because of the symmetry of the bilayer, we re- 
quire invariance with respect to the transformation (0) 



H hm = ^n(X7 2 h) 2 + SDet (h,ij) -j\7 2 h (V • m) 

+ \tm 2 + X -K x (V • m) 2 + X -K 2 (V x m) 2 (23) 



H - = -Bu 2 + -A(Vn) 2 + cVu • m 

um 2 2 



It' m 2 + \k' x (V • m) 2 + \K' 2 (V x m) 2 (24) 



The total energy therefore splits up into a contribution 
involving the average shape h and the average tilt m, 
and a decoupled contribution H involving the dilation u 

1 ^ um 

and the tilt-difference m. The term with coefficient 7 > 
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Fig. 2. Coupling between the tilt-difference m = |(mi +012) 
and the thickness gradient Vu, via the term c Vu • m. 

is responsible for the ripple phase of tilted membranes fl7| , 
|l8| . Similarly, the term with coefficient c can produce a 
"ripple" instability in which a thickness modulation occurs 
together with a tilt-difference modulation [Q. From the 
tendency of the molecules to orient perpendicular to the 
chain-water interface we expect c > (Fig. @). 



2.4 Equilibrium equations and energies 

The total elastic energy of the membrane is given by 

H = H hm +H~= [d 2 rH hm + U 2 rH~. (25) 



The equilibrium membrane configuration are those mini- 
mizing TL with respect to all possible local variations of the 
structural fields. The four corresponding Euler-Lagrange 
equations, namely 8TL/8H — 0, 8TL/8xa = 0, STL/Su = 
and STL /Sin. = 0, are explicitly 

kV 4 /! = 7 V 2 (V • m) (26) 
t m - Ki V (V • m) + K 2 V x (V x m) = -7 V (V 2 /i) 

(27) 

and 

Bu- AV 2 !i = cV • m (28) 
t'm - K[ V (V • m) + K 2 Vx(Vxm) = -cVu.(29) 

The calculation of the energy of equilibrium configurations 
can be simplified in the following way. Integrating TL ^ by 
parts yields 



H ~ = - I & 2 r [u — + m • — + TL' 
2 / \ 8u 8m 



(30) 



with 



TL' - = - q> d£ n ■ [A uVu + c u m + K\ ( V • m) m 

um 2 T 

- K' 2 (Vxm)xm] , (31) 



where the last integral is restricted to the boundary of the 
integration domain, whose normal is n. For equilibrium 
configurations, TL ^ reduces to TL 1 ^ since the first term 

of ( j30| ) vanishes. This provides a very useful simplification. 

On finds similarly that 'Hhm reduces for equilibrium 
configurations to 

Km = ^fdin- [nV 2 hVh- K,hV(V 2 h) 

- 7 (V • m) V/i + 7 h V(V • m) - 7 (V 2 /i) m 
+ ifi (V ■ m) m (V x m) x m] 



(32) 



2.5 Orders of magnitude 



For biological membranes, the bending constants n > 
and k < have relatively high values ~ 10~ 12 erg (~ 
25 IcbT) (|]. The typical value of the membrane area- 
stretching coefficient k ~ 100 erg/cm 2 [|J allows to de- 
termine the dilation modulus via B = fc/(2a) 2 , where 
a ~ 20 x 10~ 8 cm is a typical monolayer thickness. This 
yields B ~ 6 x 10 14 erg/cm 4 . Therefore B ~ n/a 4 : the 
membrane has a typical energy scale given by k and a 
typical length scale given by a. In the absence of exper- 
imental measurements, the other constants have to be 
estimated by dimensional analysis. We expect A to be 
~ n/a 2 (we recall that A is independent of the mem- 
brane tension). We therefore estimate A « 25 erg/cm 2 
Next, we shall assume roughly that tilting the molecules 
by a large angle compares energetically with compressing 
the membrane by half a monolayer thickness. This yields 
t ps t' « A. Then, we expect the characteristic length 
defined by the (Ki/t) 1 ' 2 to be of order a, which implies 
Ki k, Xa 2 w 10~ 12 erg. This value correctly compares with 
the bending constant. Finally, the K'^s are expected to be 
of the same order of magnitude as the K^s, and c is di- 
mensionally expected to compare with A. 



2.6 Remarks on the validity of the truncation of the 
energy expansion 

Strictly speaking, in all the microscopic theories of mem- 
branes j|,@,[l0, 11, 1^1, it is somewhat arbitrary to truncate 
the expansion at the lowest-order in the derivatives of the 
distortion field. Indeed, since the typical energy and length 
scales of the membrane are k and a, respectively, the dis- 
tance £ on which dilation perturbations relax is expected 
to be sa a. For small distortions B«a there is no prob- 
lem in neglecting quartic terms such as the one ~ (Vw) 4 : 
this term is (u/a) 2 times smaller than the leading term 
~(Vu) 2 . 

However the term ~ (V 2 u) 2 , which we have discarded, 
might be of the same order of magnitude as the leading 
term ~ (Vu) 2 if indeed its coefficient is ~ k. The problem 
is that all the terms ~ (V™m) 2 may be also comparable 



if their coefficients are 



However, at the mi- 



croscopic scale corresponding to a, the membrane is not 
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actually a continuum and there is not much meaning in 
considering high order derivatives of the thickness. It may 
therefore be a good approximation to keep only the lead- 
ing order term. 

In any case, we expect that the lowest-order truncation 
of this continuum description will give a correct physical 
picture, at least qualitatively, of the competing trends as- 
sociated with the various elastic variables. Note also that 
the truncation may be technically correct in the vicinity 
of a transition to a more ordered Lp or L'p phase with a 
different equilibrium thickness, where B might be signifi- 
cantly reduced and, accordingly, £ larger than a. 



3 Interactions among membrane inclusions 

Biological membranes contain a large number of inclusions 
such as integral proteins. Inclusions with a conical shape 
tend to curve the membrane since the lipids orient parallel 
to the inclusion's boundary in order to fill the volume. Be- 
cause of the interference between the resulting membrane 
distortions, such inclusions are subject to long-range in- 
teractions p9|, pQ^pj| ]. Inclusions also experience "Casimir" 
forces, which are due to the modification of the membrane 
fluctuation spectrum caused by their presence. 

The short-range interactions between inclusions arise 
from the local structural changes that the latter impose on 
the membrane [0,|,|,|l^,0|l2|. For instance, since pro- 
teins have a central hydrophobic region that spans the 
hydrophobic core of the membrane, a thickness mismatch 
between the hydrophobic region of the protein and that of 
the bilayer will result in a local membrane thickness per- 
turbation. Interferences between such perturbations yield 
membrane-mediated interactions that add up to the stan- 
dard screened-electrostatic and van der Waals interactions. 



3.1 Boundary conditions 

Let us consider a membrane inclusion such as the one de- 
picted in Fig. [| It is meant to model an integral protein 
with an arbitrary shape. For the sake of simplicity, how- 
ever, we assume revolution symmetry. We suppose that 
the hydrophobic region of the inclusion has a thickness 
2H that differs from the corresponding thickness 2a in 
the bilayer. The inclusion is also assumed to have a piece- 
wise conical shape with two angles 6\ and #2 pertaining 
to each monolayer and relative to the revolution axis. 

Let us consider an undistorted reference membrane 
above which the inclusion stands at a height ho. As previ- 
ously we denote by M 1 ' and h^ the positions of the upper 
and lower membrane interfaces with respect to their equi- 
librium positions in the reference membrane. Assuming 
a strong coupling between hydrophobic parts [p|p|,^0|,pT|, 
[l2| , we require the conditions that both monolayers inter- 
faces reach the inclusion at the separation line between its 
hydrophobic and hydrophilic regions, i.e., 



r 

h^\ ro ~ho-(H-a). 



h Q +H 



(33) 
(34) 




2 a 



Fig. 3. General boundary conditions imposed by an inclusion. 

These conditions are only approximate because the posi- 
tion where the interfaces reach the inclusion is equal to ro 
only at lowest order in the deformation variables. Another 
boundary "condition" , which is not imposed but actually 
free to adjust to equilibrium, is the angle [3 at which the 
mid-membrane shape h departs from the inclusion. Call- 
ing e. r the unit vector along r, this condition is 



Vft,| ro ~ (3e r 



(35) 



If we now require that the molecules within the mem- 
brane lie parallel to the inclusion's boundary, because of 
the space- filling constraint, we have the condition 



P (1) k 
P (2) k 



-0i e r , 
6 2 e r . 



(36) 
(37) 



Note that we have implicitly assumed that the revolution 
axis of the inclusion is normal to the reference plane (x, y), 
although in the most general situation it can be a tilted 
(this tilt will be zero by symmetry in the following). 



3.1.1 Decoupled boundary conditions 

In order to make use of the equilibrium equations previ- 
ously derived, we must transform these boundary condi- 
tions into conditions involving the variables h, u, m, and 
m. From Eqs. (|l]-||) and Eqs. (||-§|), we obtain 



h \r Q 

V/zL 



m 



ho , 
/3e r , 

{p - O) e r 



(38) 
(39) 
(40) 



and 



u\ ro ~ 1t , (41) 

m| ro ~ a e r , (42) 

!(#i + 62) is the average cone angle of the 
uq = H — a is the dilation and ao = \ (#2 — 61) 
the tilt-difference set by the inclusion. It is important to 
note that these two sets of boundary conditions are de- 
coupled in the same way as the corresponding equilibrium 
equations. 



where O 
inclusion, 
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Fig. 4. An array of inclusions and its Wigner-Seitz cell. 
3.1.2 Arrays of inclusions 

Following previous works, we shall calculate the consti- 
tutive energy of arrays of inclusions. Paradoxally, it is 
easier to approximatively calculate the energy of an ar- 
ray than to calculate exactly the interaction between two 
inclusions. Since membrane-mediated interactions are not 
pairwise additive, this is the correct procedure to investi- 
gate the stability of 2D crystalline structures. 

To capture the physics of an array of inclusions, the 
standard method is to consider a single inclusion sur- 
rounded by its Wigner-Seitz cell (i.e., the unit cell made 
by the perpendicular bisectors of the bonds connecting the 
lattice sites) ||0,[ll|,[[2|. The Wigner-Seitz cell is fur- 
ther idealized by a circle of radius approximatively half 
the inclusions separation (cf. Fig. |]), and the equilibrium 
equations are solved assuming revolution symmetry, with 
boundary conditions at r = ro and r — R. When applied 
to an hexagonal lattice of inclusions, this approximation 
is quite good, as it consists in neglecting high Fourier har- 
monics of order 6, 12, etc. In a gas of inclusions, it amounts 
to considering that the first neighbors effectively screen 
the other inclusions. 



3.2 Dilation-tilt-difference induced interactions in an 
array of inclusions 

Inclusions with arbitrary shapes will in general excite all 
of the four distortion modes considered in this work. How- 
ever, since we have seen that both the equilibrium equa- 
tions and the boundary conditions are pairwise decoupled, 
one can study separately, and simply add, the effects of the 
coupled dilation and tilt-difference modes and the effects 
of the coupled shape and tilt modes. 



3.2.1 Zero dilation-tilt-difference coupling 

We focus on the dilation (u) and tilt-difference (m) modes 
and, to start with, we neglect their coupling: 



Let us consider an array of inclusions, and assume, as 
previously discussed, a perfect revolution symmetry in the 
Wigner-Seitz cell surrounding an inclusion: 



u(r) and m = a(r) e r 



(44) 



Under these conditions, the most general solution of the 
equilibrium equations (|2^-^9|) takes the form 



u(r) 
a(r) 




A 3 Ki 



Co 



A 4 h 



6> 



(45) 
(46) 



in which the Fs and the K's are modified Bessel function 
and 



(47) 
(48) 




U\r 


= ua 




= ao 




= o, 




= 0, 



are two characteristic length comparable with the mem- 
brane thickness, except close to a Lp tilted phase where 
t! might be small, or close to the main-chain transition 
where B might be small. The constants Ai's, which are 
real and dimensionless, are determined from the bound- 
ary conditions: 



(49) 
(50) 
(51) 
(52) 



with a dot indicating derivation with respect to r. The 
quantities uq and cyq are the boundary dilation and tilt- 
difference, respectively. The last two conditions are re- 
quired by symmetry on the Wigner-Seitz circle. 

Figure |5| shows a typical solution for an isolated inclu- 
sion (R — > oo) and Fig. [6] shows a typical solution corre- 
sponding to an array of interacting inclusions. These pic- 
tures sketch the membrane structure: the solid line rep- 
resents the membrane shape, i.e., the sum of the equi- 
librium monolayer thickness a and the thickness excess 
u. The dashed curve represents the amplitude of the tilt- 
difference angle a. For the sake of clarity the distortions 
have been amplified in the following way: the boundary 
angle ao, the equilibrium monolayer thickness a, and the 
boundary thickness excess uq are all normalized to 1 . 

Assuming revolution symmetry, the general distortion 
energy ([n]) within the Wigner-Seitz cell, takes the form 



7i ~=ir\Aruii + cr 

urn L 



u a 



- K[raa + K' ia 2 ] R , 



(53) 



in which several terms vanish due to the boundary condi- 
tions. After eliminating constant terms, TL ~ reduces to 

° ' um. 



0. 



(43) H ~ = H u + H- 



-TTr (A Uq u\ ro + K[ Ce 6t\. 



(54) 
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a 



10 



15 



Fig. 5. Sketch of the membrane structure around an isolated 
inclusion (distortions are amplified, see text). The inclusion 
radius is ro = 3£ a (~ 60 A), = 2 and c = 0. 







-2 



a 



Fig. 6. Membrane structure between two inclusions in the 
array (amplified distortions). Parameters are as in Fig. |^. 



This interaction, which in principle depends on a large 
number of parameters (ro, R, Uo, oto, B, A, t' and K') has 
the following scaling property: 



n 



urn \ ' ' £ ' £ 



r ^ Q M5 



«0 



u Q ^B/V 



(55) 

(56) 
(57) 



which advantageously reduces the effective numbers of pa- 
rameters. At short inclusions separations, Ti^ diverges as 
(R — r ) _1 and H u goes to a negative constant. At large 
separations, both relax exponentially. Figure fj] shows a 
typical situation where an energy minimum appears, re- 
sulting from the superposition of a dilation-induced at- 
traction, that dominates at large distances, and a tilt- 
difference-induced repulsion, that dominates at short dis- 
tances. This situation manifests itself for large values of t; u , 
for which the dilation mode has the longest range, and for 



T~im 




5 R/Z* 6 7 

Fig. 7. Normalized interaction energy per inclusion TL -~- 
vs. the inclusions separation R. The curves correspond to 
ro = 3£ a (~ 60 A), s = 2, x = 1 and c = 0. The normal- 
ized energies Ti u and TL^* correspond to the attractive dilation 
and repulsive tilt-difference contributions, respectively. 



small values of the boundary tilt-difference ctQ, for which 
the repulsion is weak 

Let us estimate the magnitude of the interaction en- 
ergy, which is given by the normalization factor ttB ro ^ a «§ 
in (pSh. We choose for the tilt-difference coherence length 
a fixed microscopic value £ Q ~ 20 A and we let for instance 
0.1 < s < 10. This assumption is based on the fact that 
close to the main chain transition should exhibit some 
degree of pretransitional divergence. For the inclusion, we 
assume a typical protein size ro = 3£ Q (~ 60 A) and a 
thickness perturbation uo — 0.2£ Q (~ 4 A). With the esti- 
mated values of the material constants given in Sec. |2.5[ 
we obtain ttB Vq ^qUq ~ (lO/s 2 )^? 1 . 

In the energy graphs depicted in Fig. [?|, the values 
x = 1 and s = 2 correspond to an inclusion boundary tilt- 
diffcrcncc angle a = {x/s)(uo/£ a ) \J\jt 1 ~ 6°. The depth 
of the energy minimum is 3 x irB ro £ a Wo ~ 7 k^T . For such 
a well, we expect that the array of inclusion will crystallize, 
the distance between the boundaries of the particles being 
then 2(R-r ) ~ 2 X 0.8£ Q ~ 35A ||. 

If we consider the inclusions radius ro as fixed, the 
interaction potential as a function of R depends only on 
the parameters x and s, as can be seen in (^5|). We have 
plotted in Fig. @ the phase diagram, in the (x, s) plane, 
for a collection of identical inclusions. Distinction is made 
between a disordered (D) gaseous state and a crystal (K) 
phase. The criterion for the latter is the existence of a 
energy minimum with a depth larger than k^T. 



3.2.2 Nonzero dilation-tilt-difference coupling 

In order to study the effect of the dilation-tilt-difference 
coupling, we now assume 



c ^ , and £ u = £ Q = £ . 



(58) 



The latter condition is a simplification, which is reason- 
able far from any membrane phase transition (with £ of 
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Fig. 8. Phase diagram for a membrane with £ a ~ 20 A con- 
taining dilation-tilt-difference inducing inclusions with radius 
ro = 3£a(~ 60 A). The coupling c is neglected. (D) disordered 
phase. (K) crystal phase, as determined from the existence of 
an energy minimum deeper that k^T. 

the order of the membrane thickness). Let us define the 
coupling's characteristic length as 



t = 



(59) 



We assume I < £, otherwise the membrane undergoes the 
microscopic "ripple" instability already mentioned 

Under the revolution symmetry conditions 
most general solution of the equilibrium equations 
is given by the real part of 




u(r) 
a(r) 



A j Ku(e^j +A 2 I (e i ^ 
Aj K,(e^0 -A 2 Ii(e^ 



W-g,(60) 



xi, (61) 



where i = 
constants, 



-1, Ai and A 2 are two dimensionless complex 



and 



(62) 



The constants Ai and A 2 are determined from the bound- 
ary conditions (49- 52|) as previously. Figures || and fib] 
show a typical solution for an isolated inclusion (i?— >oo) 
and a typical solution for interacting inclusions, respec- 
tively. The same conventions as for Figs. || and|^ are used. 

The distortion energy within the Wigner-Seitz cell is 
given by exactly the same formula ( p4[) as previously. In- 
deed, although the second term of (|53|) does not vanish 
any longer in r = r , it is constant and can be omitted in 
the interaction. Note however that this energy cannot be 
splitted any longer into pure dilation and tilt-difference 
contributions. Again, H ~ has the following scaling prop- 
erty: 

m L^/oR^ (63) 

ttB r £ «o V £ ? / 

Depending on the values of ro/£, x and 0, the interaction 
energy is either monotonically repulsive or exhibits one or 
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Fig. 9. Sketch of the membrane structure around an isolated 
inclusion (amplified distortions). The inclusion radius is ro = 
3£(~ 60 A), x = -2, and ^ = 0.75 x tt/2. 
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Fig. 10. Membrane structure between two inclusions in the 
lattice (amplified distortions). Parameters are as in Fig. ^. 



several marked minima. Figure |ll| shows a typical situa- 
tion in which two minima appear. This phenomenon man- 
ifests itself for values of the dilation-tilt-difference cou- 
pling corresponding to <fi > 0.6 x n/2, where, because of 
the vicinity of the dilation-tilt-difference "ripple" insta- 
bility, the membrane has a tendency to develop damped 
undulations. 

The magnitude of the interaction energy is now given 
by the normalization factor ttB ro £ Uq. With typically £ ~ 
20 A, and again r = 3£(~ 60 A), u = 0.2£(~ 4A), we 
obtain, wit h the values of the material constants estimated 
in Sec. ^j], the typical energy scale ttB ro £ ~ 10 k^T. 
The boundary tilt-difference angle is then given by ao = 



x uo yBjt' 
assumed in Sec. 



2j) 



10° (for A ~ t' as consistently 



Therefore, in Fig. |ll[ the depths of 
~ 25 k^T and ~ 3 k^T, respectively. 



the two minima are 
Two distinct crystals might therefore appear: one with a 
distance between the boundaries of the particles of 2(i? — 
ro) ~ 2 x 1.2 £ ~ 50 A, the other with a much larger 
separation 2(R-r ) ~ 2 x 4.5 £ ~ 180 A [||. 
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Fig. 11. Normalized interaction energy per inclusion 7i — 
vs. the inclusions separation R. The curves correspond to ro = 
3c; (~ 60 A), x = -3, and <^ = 0.75 x tt/2. 
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Fig. 12. Phase diagram for a membrane with £ ~ 20 A con- 
taining dilation-tilt-difference inducing inclusions with radius 
ro = 3£ a (~ 60 A). (D) disordered phase. (K) crystal phase, 
(K„) region where n possible distinct crystalline phases with 
different particles separations are possible. 



If we consider the inclusions radius ro as fixed, the in- 
teraction potential as a function of R depends only on the 
parameters x and 4>, as can be seen from (|63|). Figure fl2| 
shows a phase diagram in the (x,4>) plane for a collection 
of identical inclusions. The symbol (D) indicates a disor- 
dered (D) gaseous state, (K) a crystal phase, and (K„) the 
possibility of n distinct crystalline phases with different 
separation distances. Again, the criterion for any crystal 
phase is an energy minimum depth larger than k^T . 

An interesting feature of the phase diagram of Fig. [l2] 
is the asymmetry with respect to the change x <-> — x in- 
troduced by the dilation-tilt-difference coupling: crystal 
phases are more likely to occur for x<0, i.e., for a thick- 
convex inclusion (uq > and ao < 0) or for a thin- concave 
inclusion (uq < and ao > 0). This symmetry-breaking 
follows from the sign c > of the dilation-tilt-difference 
coupling, that we have assumed, in order to favor the sit- 
uation depicted in Fig. ||. 

This asymmetry can be explained by simple arguments. 
First, let us recall that if a non-zero boundary tilt-difference 




"o<0 



Fig. 13. Sketch of the inclusions that tend to form 2D crystals 
(x < 0). The dashed lines show the monolayers of the unper- 
turbed membrane. (Left) Thick-convex inclusion (uo > and 
ao<0). (Right) Thin-concave inclusion (uo<0 and a?o>0). 

ao is present, the interaction is always repulsive at short 
distances. Indeed, the tilt-difference must go from ao to 
— ao from one inclusion to the other. Conversely, the dis- 
tortion associated with the membrane dilation is attrac- 
tive since the thickness mismatch is the same on identical 
inclusions. Two situations are therefore possible: if the di- 
lation relaxes on a longer range than the tilt-difference, 
a crystal phase can occur since there is a long-range at- 
traction followed by a short-range repulsion, whereas con- 
versely the repulsion simply dominates. We therefore have 
to understand how the coupling affects the relative range 
of the dilation and tilt-difference distortions. To simplify, 
let us rewrite schematically the interaction energy (]2J) as 

Km ~ " 2 + + ua + a 2 + £ 2 d 2 , (64) 

Let us first assume uq, ao > 0, which corresponds to x > 0. 
To relax the positive dilation uo, the membrane will set 
u < 0. The term iia being then negative, it reduces the cost 
of making a gradient of u. Therefore the u distortion will 
relax on distance that is somewhat shorter than £: the at- 
tractive dilation tail retracts (see Fig. [l4]). From the point 
of view of the tilt-difference, since u<0, the coupling iia 
makes it as if the potential was of the type (a — a m ) 2 with 
a m >0. Thus, on the distance £, the tilt-difference relaxes 
only up to a m ; it therefore needs a longer distance to reach 
zero: the repulsive tilt-difference tail expands (see Fig. [l4| ). 
Then, for x>0, a disordered phase is more favored, since 
the repulsive tail dominates at large distances. 

With still Mo > 0, let us now assume ao < 0, corre- 
sponding to x < 0. Now the term iia is positive: building 
a gradient of u is more costly and therefore the attrac- 
tive dilation tail expands (see Fig. |l5|). The tilt-difference, 
however, still experiences a potential of the type (a — a m ) 2 
with a m > 0, but it now starts from a negative value ao- On 
a distance £ it would reach the equilibrium value a m >0, 
it therefore reaches zero on a distance now shorter that 
£: the repulsive tilt-difference tail retracts (cf. Fig. p"5| ). 
Thus, for x<0, a crystal phase is more favored, since the 
attractive tail dominates at large distances (and then the 
repulsive one at short distances). 
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Fig. 14. Dilation and tilt-difference distortions around an 
isolated inclusion with x > 0. The inclusion radius is ro = 
3£(~ 60 A) and the coupling corresponds to = 0.2 x tv/2. Due 
to the latter, the tilt-difference tail expands (dashed line) and 
the dilation tail retracts (solid line), thereby favoring repulsion, 
i.e., a disordered phase. 
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Fig. 15. Same as Fig. [[4] but for x < 0. Now the tilt-difference 
tail retracts and the dilation tail expands, thereby favoring 
attraction, i.e., a crystal phase. 

3.3 Shape-tilt induced interactions in an array of 
inclusions 

We now focus on the shape (h) and tilt (m) distortion 
modes induced by the inclusions. Assuming again revolu- 
tion symmetry in the Wigner-Seitz cell, 

h = h{r) and m = 9(r) e r , (65) 

the most general solution of the equilibrium equations (]2 
|27j ) takes the form 

h = (ar 2 + b)logr+cr 2 + d+Al (qr) + BK (qr) , (66) 



with 



, L 2 a qA , , qB T ^ . . 
= -4 + — Ii(gr -^—Kx(qr). 



1 

K 



(67) 



(68) 





Fig. 16. Sketch of the membrane mid-surface shape (h) be- 
tween conical inclusions. The tilt of the lipid molecules at the 
inclusion boundary is (5 — 0. 



7 



L = 



tK 

Ki 
t 



- . C2 



e e -L 2 . 



(69) 

(70) 
(71) 



We assume L < £g, i.e., 7 2 > Kk, otherwise the mem- 
brane undergoes the ripple instability of the Ppi phase, in 
which undulations and periodic tilt distortions occur [ pT| 
0. We assume also L > 0, i.e., 7 > 0, since we expect 
the molecules to tilt in such a way as to relax the splay of 
the molecules in a curved membrane. 

To simplify, let us assume strictly K\ = k. We then 
have L = fj,£,g. Hence < /i < 1 is now the only pa- 
rameter controlling the shape-tilt coupling. The six real 
unknowns a, b, c, d, A, and B are determined, according 
to the general boundary conditions (|3^-^0|) for an inclu- 
sion with average cone angle 0, by 

(72) 
(73) 
(74) 
(75) 
(76) 
(77) 

The latter three conditions are required by symmetry on 
the Wigner-Seitz circle, at which the origin of the mem- 
brane height has been chosen. After solving this system, 
the total membrane free energy has to be minimized with 
respect to the free parameters ho and [3. Assuming revolu- 
tion symmetry, the general distortion energy ( p2| ) within 
the Wigner-Seitz cell takes the form 

Hhrn = ft n(r h + h)h — k h(r h + h — h/r) 

-j(r6 + 6)h + jh(re + 9- 6/r) 

- j(rh + h)9 + Kir + 6 2 ] R , (78) 

- ro 

where all the terms taken in r = R vanish due to the 
boundary conditions. The interaction has the following 
scaling property: 



h \r 


= h , 


h\r 


= /3, 


6\r 


= (3-0, 


h\ R 


= h Q , 


h\ R 


= 0, 


0\r 


= 0. 



Tihm 7T7 ( ^0 R 
— ILhrn I A 4 ; t ) t 



(79) 
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Fig. 17. Boundary tilt 6(ro) and boundary membrane inclina- 
tion j3 as a function of the distance R between the inclusion. 
The inclusions radius is ro = 3£e(~ 60 A) and the tilt-shape 
coupling 7 is zero. 
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Fig. 18. Same as Fig. [j^, however in the presence of a strong 
shape-tilt coupling corresponding to n = 0.9. 
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Fig. 19. Normalized interaction energy per inclusion TLhm 
vs. inclusions separation R. The inclusions radius is ro = 3 £e(— 
60 A) and the tilt-shape coupling 7 = 0. The dashed curve cor- 
responds to the case where the tilt is not allowed. 



distance between the inclusions (Fig. ^8|). Thus, except 
when the inclusions are very close to one another, the 
tilt is rapidly negligible. The reason is that the tilt set 
in order to flatten the membrane is always costly form 
the point of view of the —jW 2 h(W ■ m) coupling (when 
7 > 0). Hence the coupling does not favor an expansion 
of the tilt distortion, while in the preceding section the 
dilation-tilt-difference coupling did favor an expansion of 
the tilt-diffcrcncc for x < 0, which produced spectacular 
effects. 

Let us estimate the magnitude of the interaction en- 
ergy, which is given by the normalization factor irn 2 
in ( J79| ). With £g ~ 20 A, a typical protein size rp = 3$e(~ 
60 A) and G ~ 10°, we obtain ttk 2 ~ 2.5 k B T. Figure || 
shows the interaction energy per inclusion in the case of 
zero shape-tilt coupling. The interaction is always repul- 
sive; it diverges at small separations as (R — r )~ 1 and 
tends asymptotically towards the exact form 



The results are the following. Even in the absence of 
a shape-tilt coupling, there is a trade between the shape 
and the tilt modes, which is due to the boundary condi- 
tion ([74]). The membrane tends to develop a tilt close to 
the inclusions in order to flatten its shape. The typical so- 
lution for the membrane shape resembles that sketched in 
Fig. 16. The boundary tilt relaxes on a distance ~ iq^ 1 , 
typically of order a few £e's unless li is close to 1. The 
boundary tilt is a function of the separation R between the 
inclusions. When R 3> £0, the amplitude of the boundary 
tilt 6(ro) is negligible: the membrane curvature, which is 
small, only exerts a weak torque on the tilt. Conversely, 
when the inclusions are close to contact, the boundary tilt 
is a finite fraction of the inclusions average cone angle 0. 
For K\ — k, as we have assumed, this fraction is exactly 
1/2 for (1 = (Fig.||3). 

In the presence of a strong shape-tilt coupling (fx close 
to 1), the tilt relaxes on a distance significantly shorter 
than and the boundary tilt 9(r ), which gets somewhat 
larger at contact, actually relaxes more rapidly with the 



H h = 2ttk0 2 



R 2 - 



(80) 



which can be calculated analytically by completely ne- 
glecting the tilt. As it is apparent in Fig. |l^, the tilt re- 
laxes some of the interaction energy at short inclusions 
separations. For fi — > 1, we find that 'Hhm — * 'Hh at all 
separations. The effect of the tilt is therefore negligible. 



4 Conclusions 

We have developed an elastic model for membranes that 
describes at the same level large- and short-scale distor- 
tions of the bilayer. Strictly speaking, such a continuum 
theory at a molecular scale should not be exp ecte d to give 
more than semi-quantitative results (see Sec. 2.6). Never- 
theless our hope is that the theory captures the qualita- 
tive trends of the competitions between the different elas- 
tic variables. Using a systematic expansion in the mono- 
layers profiles and tilts, we have shown that the average 



12 



Jean-Baptiste Fournier: Microscopic Membrane Elasticity 



membrane shape (h) is coupled to the average molecular 
tilt (m), both being decoupled (at lowest order) from the 
membrane dilation (it) and the difference in the monolay- 
ers tilts (m), which are coupled together. 

We have used this model to study the contribution of 
the membrane elasticity to the short- and long-range in- 
teractions among inclusions. Because the boundary condi- 
tions at a membrane inclusion are decoupled in the same 
way as the elastic variables, the interaction energy can be 
calculated as simply the sum of a dilation-tilt-difference 
contribution (it— m) and a shape-tilt contribution (h-m). 

Membrane inclusions generally have a slightly convex 
or concave hydrophobic core of thickness different from 
that of the bilayer. Such inclusions will excite the cou- 
pled dilation-tilt-difference (it-m) mode. The thickness 
mismatch creates an energetic dilation corona around the 
inclusions and yields an attraction between like inclusions: 
no extra distortion occurs when the coronas overlap since 
the boundary dilations match. The tilt-difference, how- 
ever, yields a repulsion between like inclusions: going from 
ao to — ceo, it develops a strong gradient when the coronas 
overlap. Inclusion producing no tilt-difference aggregate, 
while inclusions producing a nonzero tilt-difference either 
repel one another or favor 2D crystals. The latter situa- 
tion arises for small tilt-differences, or when the dilation 
corona extends further than the tilt-difference corona. 

When the dilation-tilt-difference coupling is large, the 
distortions in the coronas exhibit damped oscillations. This 
effect occurs because of the vicinity of a "ripple" instabil- 
ity in which both the membrane dilation and tilt-difference 
become unstable. The inter-particle potential develops then 
several minima, which implies the possible coexistence of 
different crystals of inclusions having different lattice spac- 
ings. The latters can be significantly larger than the inclu- 
sions size. The inclusions most likely to form 2D crystals 
are those with either a long-convex or a short-concave hy- 
drophobic core, i.e., those disfavored from the point of 
view of the cVit • m coupling. This is because the gra- 
dient of u being more costly, the dilation corona extends 
(favoring "long-range" attraction) , while at the same time 
the dilation corona shrinks (making the repulsion occur 
only at smaller separations). Conversely, short-convex and 
long-concave inclusions have a dominant repulsion and 
should form disordered phases. 

Membrane inclusions generally have also a slightly con- 
ical shape. Hence they excite the coupled shape-tilt (h-m.) 
mode. In first approximation, the conical shape constrains 
the membrane to depart with a contact angle relative 
to the inclusion axis. The energy stored in the curvature 
of the membrane yields a repulsion between like inclusions 
in an array that diverges at short distances as and 
fall off as Br 2 . This is a many body effect, since the inter- 
action between a pair of inclusions falls off more rapidly, 
as R~ A |ig| . In the latter case, the inclusions axes rotate 
away from one another in order to minimize the curvature 
energy of the membrane. In an array of inclusions this 
rotation is zero by symmetry. 

If we allow for a tilt of the lipids, the membrane can 
depart with a smaller contact angle (3. In order to remain 



parallel to the inclusions boundaries, the lipids tilt then by 
(3 — 0. When the inclusions are far apart, the tilt is com- 
pletely negligible since the torque exerted by the mem- 
brane curvature on the tilt is weak. Conversely, when the 
inclusions are distant a few times the membrane thickness, 
the tilt becomes a finite fraction of 0. The interaction 
energy is then reduced, however there is no qualitative 
change in the interaction potential. As for the shape-tilt 
coupling — 7 V 2 h (V • m), it reduces the relaxation length 
of the tilt and simply reduces its effects. The reason is 
that the tilt set in order to flatten the membrane is al- 
ways costly form the point of view of the coupling, for the 
expected positive sign of 7. Hence the tilt does not propa- 
gate far away of the inclusions in the vicinity of the ripple 
instability (where both the shape and tilt modes become 
unstable). 
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